Weak vs. approximate values in quantum state determination 
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We generalize the concept of a weak value of a quantum observable to cover arbitrary real positive 
operator measures. We show that the definition is operationally meaningful in the sense that it can 
be understood within the quantum theory of sequential measurements. We then present a detailed 
analysis of the recent experiment of Lundeen et al. [lj] concerning the reconstruction of the state 
of a photon using weak measurements. We compare their method with the reconstruction method 
through informationally complete phase space measurements and show that it lacks the generality 
of the phase space method. In particular, a completely unknown state can never be reconstructed 
using the method of weak measurements. 
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I. INTRODUCTION 

The weak value of a quantum observable, as denned 
by Aharonov et al. [2( , has always remained a somewhat 
controversial concept despite the vast amount of atten- 
tion it has received over the years (see, e.g., Q and refer- 
ences therein). The purpose of this paper is to clarify this 
topic and our goal is two-fold. On one hand, we wish to 
show that the weak values can be exhaustively explained 
within the quantum theory of sequential measurements. 
For that end we briefly recall the quantum theory of mea- 
surement (Section [TTJ) and generalize the standard mea- 
surement model such that it can be applied to arbitrary 
real observables (as positive operator measures) (Section 
IIIip . We then use the theory of sequential measurements 
to define a weak value of an observable (Section IIV|) and 
compare it with its approximative values obtained from 
the relevant standard measurement. On the other hand, 
in the view of the recent paper of Lundeen et al. [l[ we 
wish to compare the role of weak and approximate mea- 
surements in quantum state determination (Section fVT). 
We analyze the method of Lundeen et al. and compare it 
with the known method of approximate sequential mea- 
surements. 

Throughout the paper we denote by T-L the Hilbert 
space associated with a physical system, and by C(H) 
and T(H) the sets of bounded and trace class operators 
acting on H. The concepts of states, observables, and 
the statistical duality they define form the rudimentary 
frame of the description of the system: a state p is (rep- 
resented as) a positive {p > 0) trace one (tr [p] = 1) 
operator, an observable is (represented as) a normal- 
ized positive operator measure (POM) E : A — > C{H) 
(defined on a er-algebra A of subsets of a set ft of 
the values of the observable), the probability measure 
A 3 X H> PppQ = tr[pE(X)] giving the measurement 



outcome statistics for the observable E in the state p. 
Usually the value space (f2, A) of an observable is just the 
real Borel space (R,23(R)), which we use subsequently. 
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For a pure state (p = p ) we use also a unit vector (vec- 
tor state) representation ip £ %, pip = ip. 



II. QUANTUM MEASUREMENTS 

Quantum theory of measurement operates on a hier- 
archy of three levels of description reflecting the options 
of restricting one's attention to the outcome probabili- 
ties at the level of the measured system, or taking into 
account system's conditional state changes, or adopting 
the most comprehensive level of modelling the interac- 
tion and information transfer between the system and 
the probe. We shall briefly describe these three levels 
since all of them play a role in our analysis, for a more 
comprehensive survey, see e.g. j4|. 

The crudest statistical level is the one described in In- 
troduction. On the next level, one also describes the 
conditional state changes of the system due to a mea- 
surement, conditioned with respect to the pointer val- 
ues. These changes are most conveniently captured in 
the concept of an instrument that is, an operation 

valued measure 1 : B(R) — > £(T(H)). This means that 
each X £ B(R) determines a positive contractive linear 
map 1(X) : T(H) -> T{U) and T{X)(p) is the unnor- 
malizcd conditional output state. Each instrument deter- 
mines uniquely the associated observable via the formula 
tr[pE(X)] = tr[l(X)(p)], or equivalently via the dual in- 
strument as E(X) =1(X)*(I). 

The most detailed descriptions of measurements are 
given by the measurement schemes where the coupling 
of the physical system and the probe (aka measuring 
apparatus) is considered. More precisely, a measure- 
ment scheme is a 4-tuple M. = (/C, <t, $, Z), where JC 
is the Hilbert space of the probe, a its initial state, 
$ : T(H ® K) — > T(H <8> K) a state transformation (i.e. 
a trace-preserving operation) modelling the interaction 
between the system and the probe, and Z : B(M) — > 
C(JC) the pointer observable (POM). Each measurement 
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scheme determines a unique instrument via 

l(X)(p)=tr c [$(p®(j)/®Z(X)] (1) 

where trjc[-] denotes partial trace over IC. As a conse- 
quence, a measurement scheme also determines uniquely 
the measured observable E, and one has 

tr[pE(X)] = ti[X(X)(p)} = tr[$(p <g> fi)/ <g> Z(X)]. (2) 

In most cases $ can be given by a unitary operator [/ on 
T-L <g> /C so that $(p (g) cr) = Z7(p (g) cr)f7*. It is often conve- 
nient to allow the pointer observable Z to have different 
values than the measured observable E has; in such a 
case one needs to introduce a (measurable) pointer func- 
tion / : ffi — > K and adjust the conditions ([T]) and ([2]) 
accordingly, that is, replacing the set X with f^ 1 (X) for 
the pointer observable Z. The fundamental result of the 
quantum measurement theory implies that each observ- 
able E admits a measurement scheme in the sense of Eq. 
© 0. 

We wish to emphasize that the correspondences be- 
tween measurement schemes and instruments, and in- 
struments and observables are many-to-one, reflecting 
the obvious fact that a given observable may be mea- 
sured in various ways, as well as that a given instrument 
may be realized by various measurement schemes. 

Now suppose that one wants to measure a pair of ob- 
servables Ei and E2 by performing their measurements se- 
quentially, for instance, first Ei and then E2. On the level 
of instruments this leads to the sequentially composed 
instrument X i2 defined by 1\2(X X Y) = 1 2 (Y) o I\{X) 

H[ . On the statistical level such a sequential measure- 
ment defines the (sequential joint) observable M given 
by U{X xY)= Ii(X)*(Z 2 (r)*(/)). In particular, the 
Cartesian margins of M are 

Mi(X) = M(X x M) = Ei(X), 
M 2 (Y) = M(R x Y) =Ji(M)*(E 2 (y)) 

which shows the characteristic quantum feature that the 
first measurement typically disturbs the subsequent one: 
the first margin is the observable measured first, the sec- 
ond margin is a smeared (or disturbed) version of the 
observable measured second, smearing depending on the 
first measurement. It is to be emphasized that the struc- 
ture of a sequential joint observable M does not depend 
on any details of the second measurement. 

III. GENERALIZED STANDARD MODEL 

As an example (and for later use) we present a gener- 
alization of the standard model of a measurement. This 
model goes back to von Neumann [9(, for a survey, see 

Let E : B(M.) — > CW) be an observable to be measured. 
By a special version Q of the Naimark dilation theorem 
there exists a Hilbert space "Hqi & unit vector ip G %a 



and a spectral measure P A : £>(R) — > C{H ®K), with the 
corresponding selfadjoint operator A, such that 

e(x) = v;p a (x)v^ 

where V^-.Ti—t W^Wo is the embedding V^(ip) = ip^tp, 

(pen. 

Let JC = L 2 (R) be the Hilbert space of the probe. For 
each A > define the state transformation $ A : T(T-L (g 
IC) T(H ® JC) via 

$ A (p® a) = tx Ho [e- iXA ® p p® IV-X^I ® ae tXmp ] 

where P is the momentum operator in K,. Since P gen- 
erates spatial translations, it is natural to choose as the 
pointer observable the position of the probe, represented 
by the spectral measure : B(K.) — >• C{K) of the posi- 
tion operator Q. Due to the coupling constant A it is now 
convenient to choose a pointer function f x (x) = \~ 1 x. 
The 5-tuple M x = (IC,a, $ A , P Q ,f x ) constitutes a mea- 
surement scheme with the intention to measure the sys- 
tem observable E. We call this a generalized standard 
model for E. 

The instrument as well as the observable actually mea- 
sured are now easily computed from (JT|) and © with the 
adoption of the pointer function f x . Since we consider 
their A-dependence we explicitly parametrize the instru- 
ment as well as the observable by that only. For nota- 
tional simplicity, we assume that the initial probe state 
cr is a pure state given by a unit vector 4>. One obtains 
the associated instrument and its dual 

I x (X)(p) = I ti no [K x V^pV;K* x ]dx, p G T(H), 
J x 

1 X {X)*{B) = { V;K* x (B®I)K x V,pdx, BeC(H), 
Jx 

where 

K x = VA<H-A(A - x)) = J VA<H~A(y - x)) dP A (y), 

for all x G M. The observable actually measured by this 
scheme is a smeared version [i x * E of the desired one, 
where the convolution is defined as 

H x * E(X) = J p x (X - x) dE(x), X G B(R), 

and the convolving probability measure is defined 
through p x (X) = (<f>\pQ(\X)<t>). 

We wish to emphasize that although the generalized 
standard measurement model as well as the associated 
instrument depend on the used Naimark dilation of E 
(aka ancilla) the actually measured observable p x * E is 
independent of the used dilation. In particular, this re- 
sult remains tacit if E is a spectral measure (in which 
case no ancilla is needed). 

The (generalized) standard model for an observable 
E constitutes its approximate measurement in the sense 
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that the actually measured observable /i A * E is a convo- 
lution of E with the probability measure /i A . The approx- 
imation depends on fi x and the degree of approximation 
can be quantified in many different ways, for instance, by 
the standard deviation of fi X . 

The set of possible values, that is, the set of pos- 
sible measurement outcomes of the observable /i A * E, 
namely the support of the POM, supp(^ A * E), is typ- 
ically much bigger than the set of possible values of E, 
supp(E). For instance, if E is a two- valued spectral mea- 
sure, supp(/i A * E) can be anything from this two-point 
set till the whole of R, depending on the support of /j, x . 
We call the numbers in supp(/it A * E) the approximate val- 
ues of E obtained from the (generalized) standard model 
Ai x . We also note that if one chooses the 'smearing' 
measure /j, x such that its expectation (or average) value 
is zero then the actually measured observable /i A * E and 
the observable E are statistically indistinguishable in the 
level of first moments (expectation values). However, the 
actual observables (POMs) coincide if and only if fi x is 
a point measure concentrated at the origin, a condition 
which is clearly impossible for our specific measurement 
scheme A4 X . 



IV. GENERALIZED WEAK VALUES 

The sequential measurement scheme can be used to 
define what Aharonov et al. called the weak value of an 
observable [2( . The intuitive idea behind the weak value 
is that by letting the strength of the interaction between 
the object and probe become sufficiently weak, the dis- 
turbance caused by the first measurement on the system 
becomes negligible. The price to be paid is that the first 
measurement becomes very poor. In other words, the ob- 
servable becomes more and more smeared. However, by 
a clever choice of the probe state it is possible to control 
the measurement so that the average value of the first 
measurement remains the same. 

The original line of reasoning can be generalized to 
cover arbitrary pairs of observables. In order to make it 
rigorous, we need some technical details. Let E : B(M) —> 
C(H) be an observable. For all ip,tp E% denote by E^ :V> 
the complex measure X h-> (ip\E(X)tp). Let T>(x, E) C H 
denote the subspace of those tp for which the identity 
map x i— > x is E^ jV -integrable for all ip E %. The first 
moment operator of E is then the linear operator E[l] : 
T>(x, E) —> H defined as 

(ip\E[l]tp) = J xdE^x), if E V{x,E),ip E H. 

The domain of E[l] contains as a subspace the square- 
integrability domain T>(x, E) consisting of those tp E H for 
which the function x i— > x 2 is integrable with respect to 
the positive measure E VtV . In general these two domains 
are different but it may happen that they coincide. This 
is the case for example for spectral measures, i.e., sharp 
observables in which case we simply have P^l] = A. 



Definition 1. Let E, F : B(R) -> C{H) be observables, 
Lp E V{x,E), \\tp\\ = I and let Y E B(R) be such that 
F(Y)tp 7^ 0. The weak value of E in a vector state tp 
conditioned by F(Y) is 



E^,F(F)) 



(y|F(Y)E[l]y) 

(<p\HY)<p) 



(3) 



Notice that if E is of the form \i * P for a probability 
measure fx, with = J xdfi(x) = 0, so that E[l] = A, 
and F is a discrete observable with F(Y) = \rj)(r}\ for some 
unit vector rj E H, we obtain the familiar expression 



E«,(¥>,F(Y)) 



(p\F(Y)E[l]y) _ ( V \Atp) 
(<p\F(Y)tp) (ti\<p) ' 



We will now proceed to showing that the general defi- 
nition is operationally meaningful. Indeed, we will show 
that the real and imaginary parts can be obtained as con- 
ditional averages in two different sequential measurement 
schemes in the limit of zero interaction strength. 

Consider the generalized standard measurement 
scheme A4 X and suppose that after this measurement, 
realizing the observable /i A * E, we perform (an exact) 
measurement of F thus obtaining the (sequential joint) 
observable M A : B(R 2 ) -> L(U). If we then postselect 
only the values [x, y) with y E Y (for a fixed Y for which 
M x (Y)ip ^ 0) and normalize the probabilities we end up 
with a conditional probability measure 



X ^ 



{<p\l x {xy{E{Y))y) _ {if\M x (X x Y)<p) 
(v\l x m*(HY)W) HM A (FV) 



We now claim that the real part of the weak value is 
obtained by taking the average of the above measure and 
then taking the limit of zero interaction strength, that is, 



Re(E„,(^F(Y))) 



lim 



(<p\l x (dx)*{E(Y)) V ) 
^|l A (R)*(F(y)V) ' 



(4) 



First notice that lim A ->o <y|X A (R)*(F(Y))tp) = 
(ip\F(Y)ip) so it is sufficient to consider only the numer- 
ator for which we define A X (X) = (tp\l x (X)*(F(Y))ip). 

Now by assuming that E T>(Q) and ip E T>(x, E) (the 
latter being equivalent to tp ® ip E T>(A) — T>(x, P A )), 
and by using the translation covariance of position we 
can show that e~ iXA ® p <p <E>ip<E>(j)E V(x, I ® I ® P Q ) for 
all A > 0. It follows that the required average value is 
given by 



A A [1 



-5<- 



tp®tp®(j)\l®I®Qe 



-i\A®P 



tp®tp®4>) 



If we then make the assumptions that ((f>\Q(f>) — and 
(f) E V{QP) n V(PQ) we can calculate the limit A — > 0. 
Using again the translation covariance of position, we 
can actually calculate the limit in a rigorous manner and 
obtain 



lim A A [I] =i{E[l]tp\F{Y)tp){<j>\PQ<t>) 

A^O 

-z(F(1>|E[1H(0|QP0). 



(5) 
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Here we have used the fact that the operator identity 
V^AV,p = E[l] holds on the square-integrability domain 

V(x,E) [Hi. 

We immediately see that a sufficient condition for the 
above limit to be proportional to the real part of the 
weak value is that (<f>\{Q, P}4>) = 0; where {•, •} denotes 
the anticommutator. A similar condition was also found 
in the earlier paper (l2| . where it was expressed as 'van- 
ishing current density'. However, we wish to obtain the 
strict equality and therefore we pose the more restrictive 
condition (cj)\QP(j)) — |. Such a condition is satisfied 

x 2 

for instance by a Gaussian 6(x] = , 1 e~ ->a 2 , Under 

these assumptions Eq. (jj) is clearly valid which shows 
that at least the real part of the (generalized) weak value 
is accessible via measurements. We summarize these con- 
siderations in the form of a proposition. 

Proposition 1. Letl x be the instrument defined by the 
measurement scheme M. x — (K,, cf>, $ A , P^, f x ) where <fi £ 
V(QP)r)V(PQ) is such that {(j)\Q<j)) = and {(t>\QP<j>) = 
^ . Then 

t {<p\l x (dxy(F(Y))<p) = (y|F(r)E[l]y>) 
aToJ X { v \T\W)*{t{Y))ip) ° (<p\f(Y)<p) 

(6) 

for all tp £ T>(x, E). 

In order to obtain the imaginary part we need another 
measurement setup. For that consider a measurement 
scheme Af x = (/C, (j), $ A , P p , / A ), which differs from M x 
by the pointer observable: instead of monitoring shifts in 
the probe's position we now observe the boosts the probe 
obtains. Note that at this point the initial probe state is 
also arbitrary. This change has a significant effect on the 
scheme. In particular, the measured observable becomes 
a trivial one, X H> (<fi\P p (\X)<p)I , meaning that no in- 
formation is gained about the object system. However, 
this does not mean that the measurement does not affect 
the state of the system. Indeed, the instrument J x asso- 
ciated with this scheme is nontrivial as can be seen from 
the dual form 

J X (X)*(B)= f \${x)\ 2 V:e iXxA {B^I)e- iXxA V^dx, 

for all B £ C{H). 

Now suppose that we again perform a sequential mea- 
surement where the second observable is F. Then the cal- 
culations for the conditional average value and the limit 
of zero interaction are performed as before. In partic- 
ular, by denoting T X (X) = (<p\ J x (X)*{F(Y))tp) and by 
assuming that <fi £ T>(P) with (cf)\P<f>) =0we find that 

lim T A [I] = 2 {(j)\P 2 4>) •Im(F(F)dE[L» 

so that by assuming that (0|P 2 </>) = \ we get the desired 
result. For the Gaussians we have (4>\P 2 (f)) — j^z so that 
the above condition is satisfied with the choice A 2 = | 



for the variance. Once again we summarize this as a 
proposition. 

Proposition 2. Let J x be the instrument defined by 
the measurement scheme Af x — (/C, <j>, $ A , P p , f x ) where 
4> £ V(P 2 ) is such that (<j)\P(j>) = and {(j)\P 2 4>) = \. 
Then 

f {<p\j x (dxy(F(Y))<p) = (y|F(y)E[l]y) 

aTo7 x (<p\j*(p)*(F(y))<p) m (<p\HY)<fi) 

(7) 

for all ip £ T>(x, E). 

We have shown that the notion of the weak value of an 
observable as given in Definition [T] is operationally mean- 
ingful in the sense that there exist sequential measure- 
ment schemes which give the real and imaginary parts 
as conditional averages in the limit of zero interaction 
strength. 

V. STATE RECONSTRUCTION METHODS 

A. Reconstructing the wavefunction via weak 
measurements 

In [l| Lundeen et al. reported an experiment where 
they claim to have measured the pointwise values of 
the wavefunction of a photon using weak measurements. 
Leaving aside the claims concerning the localization of 
a photon, the state determination method nevertheless 
deserves attention. To avoid the question of photon lo- 
calization, we consider a spin-i particle. In the original 
experiment the part of the inner degree of freedom was 
played by the polarization of the photon. 

Consider the position of the particle in, say, the z- 
direction so that the spacial part of the Hilbert space 
can be taken to be H = L 2 (R). Now as a 'probe' we 
take the spin-degree of freedom, JC = C 2 , with the initial 
state |0). In order to measure the pointwise value of the 
wavefunction we divide the position space into disjoint 
intervals (/,) with the assumption that the intervals are 
of equal length and the center is labelled by xi. For each 
interval we will then perform a weak measurement of the 
two-valued observable 1 H> Qi = pQ(Ii),0 H> I — Qi — 
P^(R \ Ii) thus scanning the whole position space. In 
order to accomplish the weak measurement, the position 
and spin of the particle are coupled via the unitary trans- 
formation 

In particular, an initial vector state if (S> |0) evolves into 
a superposition 

= Qitp® (cosa|0) +sina|l)) + (I - Qi)ip ® |0). (8) 

As the pointer observable we choose either a x or o~ y or, 
more precisely, their two- valued spectral measures P CTj , 
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j = x,y. We have thus arrived at the measurement 
schemes M) = (K, |0), $ Q , P^), i € N, j = x,y. 

After the first measurement we measure the momen- 
P p of the system and postselect the values which 



turn 



lie in the small interval J e = (- 
probabilities are then given by 



2' 2 



). The conditional 



{±1}^ 



(* Q |P p (J e ) ® P^({±l})* a ) 

(* Q |P P (J e ) ®I^ a ) ' 



(9) 



Now if we calculate the conditional average, and perform 
additional scaling by the factor 2 sin a, then in the weak 
limit a — > we arrive at the two values 



& = Rc(^|P p (J e )Q.^) 
TH = lm(^\P p {J e )Qiip) 



(10) 
(11) 



again the standard measurement scheme A4 X with the 
exception that we take the position P Q to be the observ- 
able we wish to measure so that no dilation is needed and 
the transformation <I> A becomes the usual unitary trans- 
formation. The observable actually measured is thus a 
smeared position /i A * P Q . 

Suppose that after the measurement Ai x a momentum 
measurement is performed. This results in a sequential 
joint observable M : B(M. 2 ) —> C(H) which is in fact a co- 
variant phase space observable. The structure and prop- 
erties of such observables have been studied extensively 
and are well understood [H|,[l4j]. Indeed, the observable 
M is neatly presented as 



\(Z) = — I 



W qp \4. x ){<t> x \W*dqdp 



(13) 



where £j and rji refer to the measurements of a x and 
<j y , respectively, with the fixed position interval Ii. The 
claim of [l[ is that the quantities (flOl) and (fTTj) give the 
real and imaginary parts of the wavefunction at the point 



Xi, that is, 



(p(xi) ~ constant • (ip\P p {J e )Qi<p) ■ 



(12) 



This is the more precise meaning of Eq. (2) and (3) in 
It is now straightforward to show that for a sufficiently 
regular tp, such as a compactly supported C°°-function 
with <p(0) 0, Eq. (fl2|) is true when the lengths of the 
intervals approach zero. 

The proposed method has obvious limitations as a 
method of state determination. Even after scanning 
through all the disjoint intervals Ii, the method can suc- 
ceed only for those state vectors <p for which P p ( J e )f ^ 0. 
Indeed, if the momentum of the system is localized out- 
side the vicinity of the origin then P p ( J e )<p = and no in- 
formation can be obtained from the measurement. Thus, 
this method is far from being generally valid. Moreover, 
if it happens for instance that P p {J t )<p = <p, then it is 
known from the basic Fourier theory that all the com- 
ponent vectors Qi<p, i £ N, are nonzero and one has to 
scan through all the intervals Ii in order to determine the 
state. 



B. Informationally complete sequential 
measurements 

As an alternative, and completely general method of 
state reconstruction we present the more common ap- 
proximate sequential measurement of position and mo- 
mentum where the weak limit is not taken. Consider 



where 4> x (x) = \f\cf>(—\x) and W„n = e 1 ^ e 
the Weyl operators. As proved in [15|, if the initial probe 
state satisfies {(j> x \W qp 4) X ) ^ for almost all (q,p) € M 2 , 
then the observable M is informationally complete [l6j], 
meaning that the initial state of the system is uniquely 
determined by the measurement outcome statistics of this 
observable. This implies that by choosing the probe state 
to be for instance a Gaussian, it is possible to uniquely 
determine the state of the object system with a single 
measurement scheme. A well-known optical implementa- 
tion of this measurement scheme is provided by an 8-port 
homodyne detector with a strong local oscillator. For a 
detailed analysis, see e.g. [l7| . 



VI. CONCLUSIONS 

We have given a general definition of a weak value of 
a quantum observable and shown how it is related to 
sequential measurements of pairs of observables. In par- 
ticular, we have shown that the real and imaginary parts 
of the weak value can be obtained as limits of conditional 
averages. We have also considered the determination of 
an unknown quantum state using both weak and approx- 
imate measurements. We have seen that the state recon- 
struction method of [l[ lacks the generality of the known 
methods of phase space measurements. 
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